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Qh! Abstract 

. We consider possible dynamical models for a light fermion confined by a 

potential field. With the Dirac equation only Lorentz scalar confinement yields 

- T— I I 

^ ' normalizable wavefunctions, while with the "no pair" variant of the Dirac 

Cu I equation only Lorentz vector confinement has normal Regge behaviour. A 

systematic investigation of Regge properties and phenomenological properties 
is carried out, including calculations of the Isgur-Wise function. We point out 
that the Isgur-Wise function provides a sensitive test of confinement models. 
In particular, the slope of the IW function at zero recoil point is found to be 
^'(1) ~ —0.90 for the Dirac equation with scalar confinement, and ^'(1) — 
— 1.20 for the no pair equation with vector confinement. Using heavy-light 
data alone we argue against scalar confinement. 



1 Introduction 

The description of the dynamical confinement of a hght fermion by a central force 
requires an examination of various wave equations. For example, the direct applica- 
tion of the Dirac equation with Lorentz vector confinement leads to pair production 
on hadronic time scales which is not observed. In this case the "no pair" (NP) or 
"quenched" variant of the Dirac equation is used to obtain a physical and "well 
posed" result. 

In this paper we survey the solutions of the Dirac and NP equation for both 
Lorentz vector and scalar confining potentials. We use the results in each case 
to compare to heavy-light meson data. Dirac scalar confinement and NP vector 
confinement are capable of accounting for the measured meson masses as well as the 
Isgur-Wise function. The values of the parameters resulting from the fits however 
provide interesting clues to the correct form of the confinement. In particular, we 
will argue against scalar confinement on the basis of heavy-light data alone. 

We begin in Section ^ with a brief discussion of one particle fermionic wave 
equations and their properties. Classical arguments reveal the Regge behaviour and 
categorize the effects of the negative energy states on the positive energy solutions. 
In Section ^ we discuss our numerical method and reproduce the pure Coulomb NP 
solution of Hardekopf and Sucher We then successivly consider scalar and vector 
confined mesons with Dirac and NP equations in Section § and we summarize our 
findings in Section 

2 Wave equations 

The time independent Dirac equation for a quark of mass m moving in a Lorentz 
scalar potential S{r) and the time component of a Lorentz vector potential V{r) is 




(1) 
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where the free particle Dirac Hamiltonian is 



Ho = Cfp + (3m , (2) 

and 

P(r) = V{r) + (3S{r) . (3) 

E is the energy of the "hght degrees of freedom" , and a and (3 are the usual Dirac 
matrices p|. 

When one attempts to use the Dirac equation to solve the Helium atom one finds 
[]T], ^ that no normalizable solutions are possible. This "continuum dissociation" is 
a form of the Klein paradox when the negative energy states mix with the positive 
energy states. Since Helium clearly exists, one must rephrase the Dirac equation to 
suppress this mixing. The result is the NP equation 

An analogous phenomenon occurs for a single particle moving in an increasing 
(confining) Lorentz vector potential. For a very long time it has been realized 
that there are no normalizable solutions to the Dirac equation in this case. This 
again is an example of the Klein paradox [Q. For a potential which rises sufficiently 
fast and a sufficiently small quark mass, the states corresponding to the free negative 
and positive energies mix resulting in the loss of normalizability. In both these cases 
a simple alternation to the Dirac equation avois this mixing. The free Dirac Casimir 
projection operators are 

A. = ^^, (4) 

where 

Eo = v/p2 + m2 . (5) 
These projection operators have the well known properties, 

= A± , 
A+A_ = A_A+ = , 



3 



A+ + A_ = 1 , 



(6) 



If A_|_ acts on the Dirac equation (P from the left we obtain by use of @ 

A+{Ho + P-E)^= [HoA+ + A+P(A+ + A_) - EA+]^ = . (7) 

The A-i/j term represents pair production in the interaction [|l|. The NP approxi- 
mation is obtained by dropping the A-ip term. Defining a new wave function 

= A+V , (8) 

the resuhing NP equation is 

{Ho + P+~E)(P = , (9) 

where 

P+ = A+FA+ . (10) 
2.1 Classical turning points 

The s-wave classical turning points provide valuable insight into the structure of 
these wave equations. By inspection, some vital properties of the solutions will be 
evident. The s-wave turning points are defined by p = in (P and (H). The turning 
point condition for the Dirac equation (P is 

I3{m + S) + V = E. (11) 

There are two turning points; one bounding the positive energy states (/3 = 1) and 
the other bounding the negative energy states when /3 = — 1. They are 

E+ = m + S{r) + V{r) , (12) 
E_ = -m- S{r) + V{r) . (13) 
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We consider in turn scalar confinement [y = 0, S = ar] and vector confinement 
[V = ar, S = 0]. The s-wave turning points given in (|r^) and ( [TB| ) are plotted 
in Figs. |l| and ^ for scalar and vector confinement respectively. The two cases 
are dramatically different. The scalar confinement turning points insure separation 
of the positive and negative state regions resulting in a well defined mathematical 
problem. The vector confinement turning points plotted in Fig. |^ allow the free 
negative energy states to rise into the positive energy region causing mixing for a 
given positive energy. This mixing is known as the Klein paradox |Q. In principle 
this means there are no normalizable bound states although for heavy quarks the 
mixing effect is negligible. 

At the s-wave turning point we have 

mp^o)^(;«). (14) 

and the NP turning points from (P) are given by 

E+ = m + S{r) + V{r) , (15) 
= -m . (16) 

These are plotted in Fig. ^ and are the same for scalar and vector confinement. 
From this point of view vector confinement has become a well defined concept since 
normalizable bound states are now expected to exist. 



2.2 Spherical solutions 

Because of spherical symmetry, we look for the solution in the form 



^ r| I , (17) 
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where 

r = I ° 1 , (18) 

and are the spherical spinors. Using the identity 

o- ■ pyf^ = iy^D^ , (19) 

where 

^± = ±- + (;f + -), (20) 
r dr r 

and 

D- = -Di , (21) 

which can be verified by integration by parts. We commute the Y matrix to the left 
to obtain the standard coupled radial equations, 

{ko + P- El) I M = , (22) 



where 



9' 



m -D_ 
—m 

V + S 

v^-^ 



(23) 
(24) 



The quantum number k labels the meson and is defined by 

k = ±(J + I) , (25) 

where the + sign means / = j + | or k = I, the — sign means I = i — \ or k — — (/+1), 
and hence /(/ + 1) = k{k + 1). 
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An analogous pair of coupled first order radial equations, can be obtained for 
the NP equation (01). We define 



A -^t±J!l (26) 



where E^^ is to be evaluated using orbital angular momentum / according to ([25|), 



I.e. 



E'o = pi + + . (27) 



The identity 



with 



A+r = YL , (28) 



/ 



L 

then gives the radial NP equation 



(29) 



■> 3 



jLo + LPL- El) I = . (30) 

9 J J 

It should be noted from (plf ) and (^) that = L. 

The radial equations for the Dirac case (p2D or the NP case (^) are solved in a 
similar manner as discussed in Section |^. The NP equation can of course be thought 
of as a Dirac equation with a coordinate and momentum dependent potential. 



2.3 Regge behaviour 

In the ultra-relativistic limit one aspect of the Dirac solutions can be inferred im- 
mediately without any detailed calculation. For large orbital angular momenta the 
Regge trajectories become linear for linear confinement and have a slope that is 
characteristic of the Lorentz nature of the confinement. The result that the sin- 
gle light quark Regge slopes are exactly double the corresponding two light quark 
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Regge slopes should be noted 0, where in each case the energy in the Regge slope 
is identified as the excitation energy (i.e., meson energy minus heavy quark mass). 

We consider high rotation |A;| ^ 1 and nearly circular orbits. The radial Dirac 
equation (p2D then implies 



9j_V + S-E_ ^ 

which becomes 



/f ^ V-S-E 



(31) 



{E-Vf-S^ = ^. (32) 



The state of the lowest energy for a fixed k satisfies = 0, or 



2 



{E - VyV + SS' . (33) 



r 



The two cases we are considering are 

• Scalar confinement: S = ar, V = 0, for which 



a' = \^ = — . (34) 



Vector confinement: 5 = 0, V = ar, giving a Regge slope of 



= ^ = 1 . (35) 

E^ Aa ^ ^ 



The Regge behaviour of the NP equation proceeds similarly. By pO|), (p^) and 
( P7D the high |/c| ^ 1 limit implies that 

m\»i) = l(^^^ ^ j , (36) 

where /c = and the NP equation potential term is 

LPL . if ' Mf''^^ ° 

4^\k 1 J \ V-S I \ k 1 
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We first note that the scalar potential cancels in this limit and hence the universal 
Regge behaviour is lost with NP scalar confinement. The radial equation ( pUD then 
yields 

-E + \V k(}^ + \V) 



The solutions are E = and a positive solution 



(38) 



E = V{r) + . (39) 
For the positive solution the minimum E for a fixed k condition then gives 



This is identical to the Dirac vector confinement slope 

We note that although the origin of the NP equation lies in nearly non-relativistic 
atomic physics it retains vector confinement in the ultra-relativistic limit. The status 
of scalar confinement is drastically different. The situation here is reminiscent of the 
scalar confinement catastrophe which occurs in the momentum space formulation 
M where confinement is lost. 



3 Description of the numerical method 

Instead of solving (^) by integrating differential equations, we choose a variational 
(Galerkin) method We expand fj and Qj in terms of a complete set of basis 
states 4>{r), and truncate the expansion to the lowest basis functions, i. e. 

N 

f-{r) ^ Y: ^l!^<l>nKk) , (41) 

n=l 

and 

N 

9jir) ^ Y ■ (42) 



n=l 
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Substituting these expressions into ( p2[ ) or (|30|), and then multiplying from the left 
by 

J r'^dr(j)mi{k) , (43) 
we end up with 2N x 2N matrix equation in the form 

Here, Hij are symmetric N x N matrices, and H12 = H^i- Diagonalizing the Hamil- 
tonian matrix yields energies and eigenvectors in terms of the basis states. The 
lowest N energies and eigenvectors correspond to negative energy states, while the 
+ 1 to 2N states describe positive energy states. Basis states, and therefore 
energies and eigenvectors of H, depend on the variational size parameter f3. The 
dependence on f3 of the lowest positive energy states (as well as the highest neg- 
ative energy states) should vanish with increasing number of basis states used. 
The pseudo-Coulombic basis states and all matrix elements used are described in 
Appendix |^. 

The Galerkin finite basis method will approach the true eigenvalues from above 
if the Hamiltonian matrix is positive definite. Although the Hamiltonians which we 
are considering here are not positive definite they still exhibit plateaus in (3, which 
become wider as the number of basis states increases. In those cases where analytic 
or alternative numerical solutions are available our Galerkin plateaus correspond to 
the correct eigenvalues. 

In order to see how the variational method works with the Dirac equation, we 
first solve pure Coulomb potential. 

Sir) = 0, (45) 
V{r) = ~, (46) 

for which the analytical solution is known 0] . In Fig. ^ we show dependence of the 
three lowest positive energy states (for m = 1 GeV, k = 0.5, k = —1 and j = ^) on 
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the variational parameter p. The full lines represent the exact analytic solution of 
the Dirac equation. As we increase the number of basis states, the plateau region 
of /? where the eigenvalues of H are the same as the exact energies enlarges. The 
variational scheme works well in this case. 

In Fig. H we illustrate the scaled energy of the ground state, 

E — m 



9 ' 



(47) 



as a function of the Coulomb constant k. The solid line shows the exact analytic 
Dirac result. 



E = mVl - , (48) 

and also our numerical solution which are the same to high accuracy, even for values 
of K close to one. The dashed line is the Coulomb NP equation ground state energy. 
As Hardekopf and Sucher []I| found the NP and Dirac solution are nearly the same 
for small k. Our NP Coulomb solution is consistent with that obtained by Hardekopf 
and Sucher [Q. 



4 Results 

4.1 Dirac equation with scalar confinement 

Properties of the spectrum of the Dirac equation with scalar confinement and a short 
range Coulomb interaction have already been investigated in 0. These authors also 
report results of their numerical calculations, which we have used as another check 
of our method. With = 25 basis states we were able to reproduce all of their 
numerical results for the eigenvalues to their given accuracy of four decimal places. 
As in 1^, we take 

S{r) = ar , 

V{r) = --, (49) 



11 



but here we go a little bit further in investigating the use of the Dirac equation with 
scalar confinement in the description of heavy-light mesons. First, we make sure 
that nothing in our final result depends on the value of variational prameter j3. In 
Fig. ^we show the three lowest positive energy states and the three highest negative 
energy states for = 0.3Gey, a = 0.2 GeV^, k = 0.5, k = —1 and j = |. Clearly, 
with = 15 we have a large region where the eigenvalues do not depend on p. 

Next, we perform a systematic fit to the observed spin averaged heavy-light 
meson states. We fix the light quark mass to be niu^a = 0.3 GeV, and vary all the 
other parameters of the model to best account for the experimental data. In Table 
^we show the results of this fit, with parameters 

mu,d = 0.300 GeV (fixed) , 
rus = 0.463 GeV , 
rUc = 1.301 GeV , 

nib = 4.639 GeV , (50) 
a = 0.308 GeV^ , 
K = 0.579 . 

As seen from the Table |I|, the agreement with experimental data is excellent, and 
values of parameters are all reasonable, except for the value of the tension a. 
From the universal Regge slope a' ~ 0.8 GeV~^, one expects a to be 

a = ^— - 0.2GeV^ , (51) 

and this is consistent with the value found from analyses of heavy onia spectroscopies 
[^. However, as we saw in (|5^), the Regge slope for the Dirac equation with scalar 
confinement is a' = ^, so that tension necessary to account for the spin averaged 
meson masses must be about a factor of | larger than the one expected from light- 
light spectroscopies and the Nambu slope (^TJ). Indeed, if we divide value of a from 
(|50D by |, we get 0.2, as we expected. 
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As we have already mentioned, Regge slope for the Dirac equation with scalar 
confinement is expected to be a' = ^. In order to verify this numerically, we fix a 
to be 0.2 GeV^, choose nig = 0, and then plot dependence of j with respect to 
As seen on the Fig. |^, slope of the Regge trajectories is one, as expected. 

Once the wave functions are known for a heavy light meson the Isgur-Wise 
function, describing the semi-leptonic B — > D^*^ decay distribution, can be evaluated 
rU|, 0. Using 



lu-l 

where 



= j,{2EJ-^t)) , (52) 



roo 

{A) = / dr r^R{r)A{r)R{r) , (53) 
Jo 

we find IW function that this model predicts. As shown on the Fig. ||, the agreement 



with ARGUS [T^ and CLEO |jT3[ data is reasonable. To calculate the slope, we use 
expression 
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e{l) = -{^ + ^E'g{r')), (54) 
For the range of light quark masses from to 350 MeV, we obtain 

^'(1) = -0.90 ±0.02 . (55) 

4.2 Dirac equation with vector confinement 

If a fermion is confined by a Lorentz vector interaction, i. e. 

S{r) = 0, 

V{r) = + ar , (56) 

the Dirac equation has no normalizable solutions. As we mentioned earlier this has 
been known to be the case for over sixty years [^]. The origin of the problem is in the 
mixing between positive and negative energy states as shown in Fig. ^ Since the 
basis wavefunctions are all normalizable in our variational method it is of interest 
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to see how this problem is manefest. In Fig. ^ we show the /3-plot with vector 
confinement for the three lowest positive and the three highest negative energy 
states. We see how the three states mix in the region of (3 where there should be a 
plateau. It is interesting to note however that if one extends the apparent plateaus 
through the level repulsions a consistent result is obtained which is not too different 
from the NP result below. The largest difference is in the ground state where Dirac 
and NP are about 50 MeV apart. 

4.3 NP equation with scalar confinement 

We observed in ( ^Tf ) that for large orbital excitation the NP scalar interaction cancels 
from the NP equation. The consequent loss of linear Regge trajectories eliminates 
any conventional discussion of meson states in terms of scalar confinement in the 
NP framework. 

4.4 NP equation with vector confinement 

There is reason to hope that the NP equation with vector confinement (|56|) will 
eliminate this mixing and hence reestablish a one particle wave equation. Comparing 
the s-wave turning point structure of the NP equation in Fig. |] with the Dirac 
equation shown in Fig. ^ we observe that the NP negative energy states should not 
seriously mix with the positive energy states. 

We proceed to the numerical solution of the NP equation ( |30D by the Galerkin 
method. In Fig. we show the dependence on /3 of the three lowest positive and 
three highest negative energy solutions with nig = 0.3 GeV, a = 0.2 GeV^, k = 
0.5, k = —1 and j = |. Again we observe a robust plateau structure which widens 
as the number of basis functions increases. 

The Regge slope of the NP equation with vector confinement was shown in (^) 
to be a' = ^. Our numerical solution agrees as shown in Fig. |TT[ In this figure 
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we illustrate the leading trajectories for the two light degrees of freedom states 
corresponding to A; = ±{j + i) for rriu^d = 0.3 GeV, a = 0.2 GeV^ and k = 0.5. We 
also show several daughter trajectories corresponding to radial excitations. 

In order to recover the universal Regge slope, the tension in this case must be 
about J times the one from equation (^1]). Therefore, we fix a to be 0.16 GeV^, 
choose rriu^ = 0.3 GeV and fit to the spin averaged heavy-light meson states. Result 
is shown in Table ^, and parameters of the fit are 

mu,d = 0.300 GeV (fixed) , 
m, = 0.600 GeV , 
rric = 1.342 GeV , 

rub = 4.679 GeV , (57) 
a = 0.157 GeV'^ (fixed) , 
K = 0.676 . 

The agreement of the fitted levels to experiment is very good and comparable to the 
Dirac scalar confinement fit of Table 

Finally we use the wavefunctions with parameters of ( p7|) to evaluate the IW 
function using (^2]). The result is shown on Fig. |1^. The slope of the Isgur-Wise 



function at the zero recoil point is evaluated using ( ^4] ) to be 

^'(1) = -1.20 ±0.03 , (58) 

where the error is estimated from the variation of the light quark mass 0.25 < rriu^d < 
0.35 GeV. We observe that this slope is significantly more negative than the one 
found from the Dirac equation with scalar confinement (0), even though the same 
set of spin averaged heavy-light meson masses were used in the fit. However, if we 
compare Figs. |^ and |12|, we see that Isgur-Wise function obtained from the NP 
equation agrees with the data a bit better than the one calculated from the Dirac 
equation with scalar confinement. 
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5 Conclusions 



We have considered here the motion of a fermion in a central field. The interaction 
that we have emphasized is scalar or vector linear confinement. We are particularly 
interested in this problem because of its application to the description of heavy- 
light mesons. The wave equations considered are the Dirac equation and the "no 
pair" equation. Although the properties of the Dirac equation are well known we 
reconsider them in the light of a confining interaction and also to serve as benchmark 
for the related NP equation. Our main results are: 

1. Dirac equation with scalar linear confinement 

This is the most straightforward confinement model for a fermion. We find 
that it is a well posed problem with a unique solution and that by adjusting 
quark masses, Coulomb constant k, and tension a, an excellent fit to heavy- 
light masses can be found. The tension found (a ~ 0.31 GeV) is larger than 
normally obtained from heavy onia fits and this may be viewed as evidence 
against scalar confinement. The origin of this discrepancy is the scalar con- 
finement Regge slope = ^- The slope is fixed by the p-wave heavy-light 
states and yields a tension about 50% larger than the usual value. If we force 
the tension to be 0.2 GeV'^ then although the fit to the energy levels is still 
good, the Isgur-Wise function becomes shallow with a slope ^'(1) ~ —0.70 and 
does not fit the experimental data well. We conclude that although scalar con- 
finement gives mathematically consistent solutions, it does not seem to agree 
well with experiment. 

2. Dirac equation with vector confinement 

There are no normalizable bound state solutions in this case. 

3. NP equation with scalar confinement 
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Scalar confinement in the no pair equation does not yield normal quasi-linear 
Regge trajectories and hence does not correspond to universal Regge be- 
haviour. 



4. NP equation with vector confinement 

This model of fermionic confinement also is viable. Vector NP confinement, 
along with an attractive Coulombic short range interaction, again gives good 
fits to the data. Vector confinement yields a Regge slope of a' = ^, one half 
of the scalar value. 

There are thus two alternative confinement models yielding linear Regge tra- 
jectories: scalar Dirac and vector NP. They both account for the data well. The 
differences though are interesting. The parameters, most notably the tension, are 
different when fitted to the data. 



As we observed in Section |2.3| the Regge slope for scalar and vector confinement 
differ by a factor of two. Since the tension appears both in s-wave dynamics and 
in rotational states experimental data will ultimately decide the correct result. The 
tension prefered by heavy onia actually lies between the scalar Dirac and vector 



NP heavy-light values. As we have pointed out |]T4|, |1^ the heavy-light Regge slope 
a' = — obtained in the relativistic flux tube model is consistent with the heavy 
onia value. 

Another piece of information which should soon shed additional light on the 
proper model of conflnement is the Isgur-Wise function. Both by experiment and 
from the lattice simulation of QCD, accurate values of the IW function (or its slope at 
zero recoil point) will be available. The Isgur-Wise function appears to depend fairly 
sensitively on the conflnement model. As we have seen, Dirac scalar conflnement 
yields = —0.9, while NP vector conflnement gives ^{1) = —1.2. Both of these 
models provide excellent flts to the same data set. 
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Finally, we should mention that we have been led into these questions by our 
investigation of the relativistic flux tube model |T^, |T3|, |T^. For low orbital angular 
momentum states the flux tube model is similar to vector confinement. It appears 
from the results presented here that a NP type equation will be appropriate for the 
flux tube with one light fermion. 



A Basis states 



The pseudo-Coulombic radial basis states are those used in previous calculations 

Ruir) = NuPH2Prye-^^Lf^'i2(3r) , (59) 

where 

In these equations we assume < z < — 1. For computational precision and 
efficiency, the matrix representation of all operators have been calculated analytically 



171 . Forp^ 



and ^ we have 



X M n^-^C^^ + + + 1)(2^ + 2/ + 3) 



(2/+ l)(2/ + 3) 



r 



1 



2i + 2l + 3)6^J 



j{j + 21 + 2)5i,,_i - Ji{t + 21 + 2)6^J+l 



l + lNu ' 



2(3' 



r'"' {l + l){2l + l){2l + 3)Nu 



[(2/ + l)(j-z) + 2j + 2/ + 3] 



(61) 

,(62) 
(63) 
(64) 
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where we assume i < j, while results for i > j can be obtained by simple reflection 
due to the symmetry of operators. The only other matrix element used was 



^dr''' l + l 



-1 



I < J 

i = j 
i > j 



(65) 
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TABLES 



Table 1: Heavy-light spin averaged states. Theoretical results are obtained from 
the Dirac equation with scalar confinement. Spin-averaged masses are calculated in 
the usual way, by taking | (|) of the triplet and j (|) of the singlet mass for the 
s(p)-waves). An estimate for the B* was taken to be 5421 MeV, in order to make 
sphttings B* — B and S* — Bg the same. 



state 


spectroscopic label spin-averaged 


q- 


n. 


theory 


error 






2^+1 L J mass (MeV) 


j 


A; 


(MeV) 


(MeV) 


CM, cd quarks 
















D (18o7) 
D* (2010) 


0" 
1- 


'So 
'S, 


IS (1974) 


1 

2 


-1 


1975 


1 


Di (2423) 
Dl (2457) 


1+ 
2+ 


'Pi ' 

'P2_ 


IP (2444) 


3 
2 


-2 


2444 





cs quarks 
















D, (1969) 
Dl (2110) 


0- 

1- 


'So' 
'S, 


IS (2075) 


1 
2 


-1 
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Table 2: Heavy-light spin averaged states. Theoretical results are obtained from 
the no pair equation with vector confinement. Spin-averaged masses are calculated 
in the usual way, by taking | (|) of the triplet and j (|) of the singlet mass for the 
s(p)-waves). An estimate for the B* was taken to be 5421 MeV, in order to make 
splittings B* — B and 5* — Bg the same. 
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FIGURES 



Figure 1: S-wave classical turning points for the Dirac equation with scalar confine- 
ment. 

Figure 2: S-wave classical turning points for the Dirac equation with vector confine- 
ment. 

Figure 3: S-wave classical turning points for the no pair equation. Note that there 
is no difference between the scalar and the vector confinement. 

Figure 4: Dependence on [3 of the three lowest positive energy states for the pure 
Coulomb case with m = 1 GeV, k = 0.5, k = —1 and J = |- Exact solutions of the 
Dirac equation are shown by full lines, while our variational solutions are shown by 
dashed lines (shorter for N — 15 and longer for N — 25 basis states used). 

Figure 5: Scaled eigenvalue e for the ground state (j = ^, k = —1). We used 
m — 1 GeV and N = 25 basis states. The no pair equation result is shown by the 
dashed line, while the Dirac equation result is shown by the full line. 

Figure 6: Dependence on /? of the three lowest positive energy states and the three 
highest negative energy states for the Dirac equation with scalar confinement, with 
ruq = 0.3 GeV, a = 0.2 GeV"^, k = 0.5, /c = — 1 and j = |. Full lines correspond to 
N — 25, and dashed lines correspond to A'' = 15 basis states used. 

Figure 7: Regge trajectories for the Dirac equation with scalar confinement. We have 
chosen rriu^d = 0, a = 0.2GeV'^, and n = 0.5. Full lines correspond to k — —{j + |), 
and dashed lines to A; = j + \- To ensure that all calculated energies are correct, we 
used N — 100 basis states, and kept first 15 states. 
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Figure 8: IW function for B decays calculated from the Dirac equation with scalar 
confinement. Values for the light quark mass rriu^d, tension a and short range poten- 
tial constant k are taken from the (|50|). For the sake of clarity, error bars are shown 
only for the CLEO data. 

Figure 9: Dependence on f3 of the three lowest positive energy states and the three 
highest negative energy states for the Dirac equation with vector confinement, with 
rUg = 0.3 GeV, a = 0.2 GeV^, k = 0.5, A; = -1 and j = |. We used = 15 basis 
states. The mixing between positive and negative states is evident. 

Figure 10: Dependence on (3 of the three lowest positive energy states and the three 
highest negative energy states for the no pair equation with vector confinement, with 
rrig = 0.3 GeV, a = 0.2 GeV"^, k = 0.5, k = —1 and j = |. Full lines correspond to 
= 25, and dashed lines correspond to = 15 basis states used. 

Figure 11: Regge trajectories for the no pair equation with vector confinement. We 
have chosen m^j ^ = 0.3 GeV, a = 0.2GeV^, and k = 0.5. Full lines correspond to 
/c = — (j + |), and dashed lines to A; = j + |. To ensure that all calculated energies 
are correct, we used iV = 50 basis states, and kept first 10 states. 

Figure 12: IW function for B decays calculated from the no pair equation with 
vector confinement. Values for the light quark mass niu^d, tension a and short range 



potential constant k, are taken from (|57|) . For the sake of clarity, error bars are 
shown only for the CLEO data. 
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